Theory of Electromotive Force Induced by Domain Wall Motion 
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We formulate a theory on the dynamics of conduction electrons in the presence of moving magnetic 
textures in ferromagnetic materials. We show that the variation of local magnetization in both space 
and time gives rise to topological fields, which induce electromotive forces on the electrons. Universal 
results are obtained for the emf induced by both transverse and vortex domain walls traveling in 
a magnetic film strip, and their measurement may provide clear characterization on the motion of 
such walls. 
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The interplay between electron transport and magnetic 
dynamics is a central problem of spintronics research. 
It has been known that the presence of a domain wall 
(DW) can change the electrical resistance of a ferromag- 
netic conductor [H, 0, IE 0] • It has also been demonstrated 
that an electric current can drive a DW through coupling 
between the conduction electrons and the local magnetic 
moments @, @, 0, IE IE EE EH EE EE EL El- The re- 
verse of this effect, i.e., electron transport induced by 
a moving DW, has been proposed by Berger [l6j] based 
on phenomenological arguments about twenty years ago. 
Recently there is renewed interest in this effect and the 
result for ID transverse domain wall has been derived 
rigorously using various approaches 17 , IE HI > connect- 
ing to Berry phase effects on electron spins in a magnetic 
texture [20, |2l|, [22| • However, realistic DW profiles and 
their motion are far more complicated than the simple 
ID model 23, E|, typically involving vortices. So a gen- 
eral microscopic theory applicable to higher dimensions 
is very much desired. 

In this Letter, we provide such a theory within the 
framework of semiclassical dynamics of electrons in a 
magnetic background which varies slowly in both space 
and time. Indeed, the width of a typical DW in a ferro- 
magnetic nanowire is about a few hundred nanometers, 
which is much larger than the electron Fermi wavelength, 
and the DW speed is much smaller than the electron 
speed. One can therefore consider the semiclassical for- 
malism with adiabatic approximation where the spin of 
the conduction electron follows the direction of local spin 
vector. We find that Berry phase terms [2(| enter into the 
equation of motion as a pair of topological fields, acting 
locally like electric and magnetic fields on the electrons. 

While our theory naturally reproduces previous results 
on ID transverse walls, we consider carefully the effect 
of a moving vortex wall. We predict both longitudinal 
and transverse voltage with universal results which may 
be used to provide clear characterization of the wall mo- 
tion. Extra voltage due to nonadiabatic effects is also 
estimated at the end of the paper, and is found to be 
subdominant. 

To construct our theory, we consider a ferromagnetic 



thin film, which is taken to be the x-y plane. The time 
evolution of local magnetization can be driven by a uni- 
form applied magnetic field (x-direction). Our model 
Hamiltonian then takes the following form, 



H = H Q [q + (e/H)A(r)] - Jn(r, t) ■ a - ha x 



(1) 



The first term is the bare Hamiltonian for a conduction 
electron, q is the Bloch wave- vector, and A(r) is the vec- 
tor potential of the external magnetic field. The second 
term is the s-d coupling between a conduction electron 
and the local d-electron spin along direction fi(r,t), and 
J is the s-d coupling strength. The last term represents 
the Zeeman coupling between electron and the external 
magnetic field, with h = \g s \iB^- 

To apply the semiclassical wave-packet formalism for 
the conduction electrons [2^], we first write down the 
local Hamiltonian at the center position of the electron 
wave-packet r c , 

H c = H [q + (e/h)A{r c )} - K{r Cl t)n'(r c , t) ■ a, (2) 

where n'(r c ,t) is the unit vector of the exchange plus 
Zeeman field, while K(r r , t) is its strength. As discussed 
by Sundaram and Niu [25|, we introduce the gauge in- 
variant crystal momentum k = q + (e/K)A.(r). For now, 
we only consider the majority carriers whose spins are 
polarized along n'(r C) i), and spin minority carriers are 
considered in the end. The position and time dependence 
of the spinor wave function gives rise to Berry curvatures 
in space and time, which can affect the dynamics of an 
electron wavepacket. More specifically, we find that the 
equations of motion for the wavepacket center are (sub- 
script c is dropped here) 
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(3) 



(4) 



where £o is the Bloch band energy obtained from Ho. 
Originated from Berry curvatures in real space and time, 
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FIG. 1: (color online). Schematic picture of a transverse do- 
main wall. The local spin direction changes within a wall 
plane which makes an angle ip with the vertical direction. 



two new fields C and D appear in the equations of mo- 
tion. They are entirely due to the spatial and temporal 
variation of local spin textures. In terms of the spherical 
angles (6, 4>) specifying the direction of n' , the fields C 
and D are given by 



C{r,t) ee i sin (V9 x V</>) , 



D(r,t) ee - sin 
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(5) 



(6) 



Field C is similar to the gyro vector used in the discussion 
of Bloch line dynamics [26(. From their appearance in 
equation ^ , we observe that C acts like a magnetic field 
while D behaves like an electric field. We note that our 
theory applies locall y to the dynamics of electrons, while 
recent work of Ref. 17, [HI focuses more on the global 
aspect of Berry phase effects. 

Next we apply these formulas to the one-dimensional 
domain wall case where and <j) only depend on co- 
ordinate x of the 2D plane. We then have C = 
and D = De x . For a transverse domain wall in the 
absence of external magnetic field, the local spins ro- 
tate within the 2D plane with the profile = ir/2 and 
4>(x) = cos -1 tanh[(x — X)/X], where X denotes the cen- 
ter of the DW and A is the DW width [2fJ. When an ex- 
ternal B field is applied, not only will the DW propagate, 
but the local spins will also get tilted out of the x-y plane. 
As shown in FigJU these spins still lie within a "wall 
plane" which makes an angle ip with the vertical direc- 
tion. Below Walker's breakdown field, this plane is fixed, 
and the angles and <j) are functions of [x — vt), which 
makes D = 0. Above Walker's breakdown field, this 
plane changes with time (with rate denoted by dip/dt) 
[26j |. which leads to a non-zero D field. 

If the system is bounded electrically, in a steady state, 
the DW motion induced adiabatic force h(dK/hdx — D) 
must be balanced by the gradient of electrochemical po- 
tential. The emf along the direction of DW motion, which 
is the measured voltage change, is given by 



V = - [ (— -D\dx = ~ (— + — 
x e / V hdx J e \ dt h 



(7) 



In the calculation, we use the condition \h\ <C \J\, which 
is usually the case in experiments on ferromagnetic ma- 
terials. The first term on the right hand side repre- 
sents the so-called AC ferro-Josephson effect proposed 
by Berger in 1986 based on phenomenological consider- 
ations [l6|. More careful analyses on this effect have 
been done recently using different approaches 17, 3, ljj ]. 
The additional term proportional to the field is due to 
the difference in Zeeman energy on the two sides of the 
DW, which should appear when we suddenly turn on the 
external field. However, since the spin relaxation time 
t s ~ 10 _12 s is much shorter than the characteristic time 
for DW motion, the voltage associated with this term 
cannot be measured in experiment 27 1. 

Depending on the thickness and the width of the sys- 
tem, a stable DW in a thin film can also take a vortex 
structure 23|, |24|. In fact, most of the experiments done 



so far involve vortex walls rather than transverse walls. 
In the following, we apply the formula developed above 
to study the case of a single vortex in a nanowire. As- 
sume that the vortex profile is characterized by a core 
radius a, which is about a few nanometers, and an outer 
radius R, which is comparable to the wire width w. For 
a vortex centered at X(t), we may approximately write 
9 = f [1 -pexp(-\r -X(t)\/a)] for \r - X\ < R and 
= 5 beyond the outer radius, p = ±1 is the polariza- 
tion. In both cases, we have <j> = garg(r — X(t)) + 
for \r — X\ < R, where q = ±1 is the vorticity of the 
vortex (g_= —1 is also referred to as antivortex in lit- 
erature [301] ). and c = ±1 indicates its chirality. For a 
steady state motion of the vortex, and <f> are functions 
of (r — vt), where v = X. Then 

D = — ~ sm0[(v ■ V(f>)\70 - (v ■ V0)V<£] = C x v, (8) 

which resembles the relation between E and B fields. 
Because a <C R ~ w, these fields are concentrated within 
the core region, where 



C = M ^-e-^co S (|e-^) 



(9) 



and D is obtained by Eq.®. An important property of 
the C field is that its total flux is a constant, i.e., 



Cd\ = I 



sin 



pqire z 



(10) 



which is topologically invariant — independent of the de- 
tailed profile of the vortex. 

The two-dimensional character of the vortex domain 
wall makes the calculation of the induced voltage a bit 
complicated. Unlike the one-dimensional transverse wall 
case, the force field D now has a curl. The gradient of 
the electrochemical potential can only cancel the longitu- 
dinal part of this force field. Therefore, we need to solve 
the Poisson equation V 2 V = (fi/e)V • D with Neumann 
boundary condition (no current leaving the sample). It 
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FIG. 2: (color online). Schematic picture of the D field for a 
moving vortex. D field is perpendicular to the vortex velocity 
and concentrated within the vortex core. 



is seen that the effect of the D field may be regarded as 
that of an electric dipole P, whose spatial extension is 
of the size of the core. The net dipole moment is equiv- 
alently given by the integral of the D field times Heo/e. 
With the help of Eq.® and ([TP]). P = pqeoir(h/e)e z x v, 
which is also a topological property of the vortex. 

Another complication arises from the Magnus force 
on a moving vortex which pushes it in the direction 
perpendicular to its velocity (28|. At low fields, this 
Magnus force is balanced by the confining potential of 
the nanowire such that steady motion is still along x- 
direction. In this case, the source term of the Poisson 
equation resembles an electric dipole pointing along the 
y-dircction, so the longitudinal voltage is expected to 
vanish in this case (baring nonadiabatic effects). 

Above the breakdown field, the confining potential can 
no longer balance the Ma gnu s force and the vortex will 
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begin a transverse motion 1J 
hits one edge, the vortex domain wall transforms into a 
transverse wall. There is a range of external field under 
which this transverse wall propagates and generates a 
voltage according to Eq. ([7} • But more probably, another 
vortex with reversed polarization will be emitted from 
the edge, travel across the wire and hit the other edge. 
These transformations continue periodically as the DW 



moves along the wire [31 



When the vortex begins transverse motion, from the 
relation D = C x v, we observe that the dipole source 
gets rotated to acquire a finite x component. This 
makes the longitudinal voltage nonzero (FiglJ). Analyt- 
ical expression for the voltage can be obtained within a 
point-dipole approximation and using the image charge 
method. The longitudinal voltage drop along the direc- 
tion of the DW motion is obtained as, 



V x = 7T , 

e w 



(11) 



with v y being the magnitude of the transverse velocity. 

This longitudinal voltage drop is proportional to the 
transverse speed, and inversely proportional to the wire 
width. The result is universal in the sense that it is hide- 
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FIG. 3: Numerical results for longitudinal DC voltage asso- 
ciated with the vortex motion. The triangular data points 
are obtained by directly solving the Poisson equation with 
v y = *yHw/ir for the sample used in Beach et al.'s experiment 
[33 ]. For comparison, a solid line is drawn based on Eq. (|12|l . 
Here also shows the mobility curve measured experimentally 
in Ref.f3l. 



pendent of the detailed wall profile including its polariza- 
tion and vorticity. In fact, as confirmed by our numerical 
calculations, this result is exact beyond the point-dipole 
approximation as long as the core is contained within 
the width and is far from the two ends of the sample. 
Measurement of this universal result may provide clear 
characterization of the vortex motion. 

The average value of this voltage depends on the aver- 
age frequency of the wall transformation. Recent exper- 
iments and simulations [3ll . [32j suggest that for narrow 
nanowires, this frequency is approximately the Larmor 
frequency. So v y /w = jH/tt, where v y is time averaged 
transverse speed, 7 is the gyromagnetic ratio and H is 
the applied magnetic field strength. In this case, 



Vr 



(12) 



which means the DC part of the voltage signal is universal 
and only depends on the applied magnetic field. This 
result is in good agreement with the voltage obtained by 
solving Poisson equation directly (Fig[3]). Here we also 
expect there are small oscillations around this DC signal 
with the Larmor frequency. 

The motion of the vortex also induces a transverse volt- 
age, which can also be calculated from the Poisson equa- 
tion. For a vortex at the center of the wire, the transverse 
voltage at the longitudinal position of the vortex is found 
to be 



Vu 



npq- 



H v x 



(13) 



which can be measured through a pair of lateral leads. 
At low fields, this voltage is a constant, and one should 
observe a pulse of transverse voltage with the above peak 
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value with its width determined by the wall speed and 
lead width. Strong oscillations occur within this pulse 
above the breakdown field when the vortex or antivortex 
executes complicated motion. 

So far we have only considered the spin majority car- 
riers and the sample being disorder free. In real situa- 
tion, for adiabatic approximation to be valid, we need 
J/(fi/r) 3> 1, where r is the carrier mean free time be- 
tween elastic scattering events [34j]. This condition en- 
sures the probability of spin flip due to collision broad- 
ening is negligible. For Permalloy thin films at room 
temperature, this condition still holds for spin majority 
carriers but breaks down for minority carriers. There- 
fore the above Berry phase effects disappear for minority 
carriers because their phases are randomized due to scat- 
tering. 

Finally, we give an estimation of the extra voltage 
due to nonadiabatic or dissipative effects. Originally ob- 
tained from force balance considerations this contri- 
bution has been related to the nonadiabatic spin transfer 
torque recently [Hl |. This voltage drop can be written 
as V™ a — 2M s R [i~ 1 v x , where Rq is the ordinary Hall 
coefficient and /i^, called the intrinsic wall mobility, is 
a measure of the electrons' contribution to the viscous 
damping force on the domain wall. For Permalloy thin 
films, M s = 8 x I0 5 A/m, R = -1.4 x l(T 10 m 3 /C and 
Hi ~ 2m 2 /C. Then V™ a is at least 10 times smaller than 
the adiabatic voltage above breakdown. However, it is 
the dominant contribution to the longitudinal voltage be- 
low breakdown. 

In summary, we have proposed a general theory for 
studying electron dynamics in the presence of moving lo- 
cal spin textures. We find that the variation of local spin 
textures gives rise to two topological fields acting on the 
conduction electrons as driving forces. Using this formal- 
ism, we reproduce the result for transverse wall motion. 
Moreover, universal results are obtained for the voltage 
induced by a moving vortex wall and its measurement can 
be used for detecting the domain wall motion. Finally, 
we estimate the nonadiabatic contributions to the voltage 
drop which shall be important below Walker breakdown. 
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